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Abstract 

These lecture notes concern recent developments in our under- 
standing of the low energy properties of QCD. Significant progress 
has been made on the lattice and the beautiful experimental results 
on the K e 4 and decays, as well as those on pionic atoms also con- 
firm the results obtained on the basis of Chiral Perturbation Theory. 
There is an exception: one of the precision experiments on decay 
is in flat contradiction with the Callan-Treiman relation. If confirmed, 
this would indicate physics beyond the Standard Model: right-handed 
quark couplings of the VP-boson, for instance. Furthermore, I discuss 
two examples where the estimates of the effective coupling constants 
based on saturation by resonances appear to fail. 

In the second part, the progress made in extending the range of 
validity of the effective theory with dispersive methods is reviewed. In 
particular, I draw attention to an exact formula, which expresses the 
mass and width of a resonance in terms of observable quantities. The 
formula removes the ambiguities inherent in the analytic continuation 
from the real axis into the complex plane, which plagued previous 
determinations of the pole positions associated with broad resonances. 
In particular, it can now be demonstrated that the lowest resonance 
of QCD carries the quantum numbers of the vacuum. 
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1 Introduction 

QCD with massless quarks is the ideal of a theory: it does not contain a 
single dimensionless free parameter. At high energies, the degrees of freedom 
occurring in the Lagrangian are suitable for a description of the phenomena, 
because the interaction among these degrees of freedom can be treated as 
a perturbation. At low energies, on the other hand, QCD reveals a rich 
spectrum of hadrons, the understanding of which is beyond the reach of 
perturbation theory. In my opinion, one of the main challenges within the 
Standard Model is to understand how an intrinsically simple beauty like QCD 
can give rise to the amazing structures observed at low energy. 

The progress achieved in understanding the low energy properties of QCD 
has been very slow. A large fraction of the papers written in this field does 
not concern QCD as such, but models that resemble it in one way or the 
other: constituent quarks, NJL-model, linear a model, hidden local symme- 
try, AdS / CFT and many others. Some of these may be viewed as simplified 
versions of QCD that do catch some of the salient features of the theory at 
the semi-quantitative level, but none provides a basis for a coherent approx- 
imation scheme that would allow us, in principle, to solve QCD. 

These lectures concern the model independent approach to the problem 
based on effective field theory, lattice methods and dispersion theory. The 
effective theory relevant for low energy QCD is referred to as Chiral Pertur- 
bation Theory, . For recent reviews of this framework in the mesonic 
sector, I refer to [l]-[5]. An update of the experimental information con- 
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cerning the effective coupling constants was given by Bijnens at the lattice 
conference in 2007 [6]. The rapidly growing information obtained in the 
framework of the lattice approach is reviewed in the report of Necco at the 
same meeting [7]. An up-to-date account of the progress made in applying 
%PT to the baryons is given in [8]. 

At low energies, the main characteristic of QCD is that the energy gap is 
remarkably small, M n ~ 140 MeV. More than 10 years before the discovery 
of QCD, Nambu [9] found out why that is so: the gap is small because the 
strong interactions have an approximate chiral symmetry. Indeed, QCD does 
have this property: for yet unknown reasons, two of the quarks happen to 
be very light. The symmetry is not perfect, but nearly so: m u and are 
tiny. The mass gap is small because the symmetry is "hidden" or "sponta- 
neously broken" : for dynamical reasons, the ground state of the theory is not 
invariant under chiral rotations, not even approximately. The spontaneous 
breakdown of an exact Lie group symmetry gives rise to strictly massless 
particles, "Goldstone bosons" . In QCD, the pions play this role: they would 
be strictly massless if m u and were zero, because the symmetry would 
then be exact. The only term in the Lagrangian of QCD that is not invariant 
under the group SU(2) xSU(2) of chiral rotations is the mass term of the two 
lightest quarks, m u uu + rriddd. This term equips the pions with a mass. 
Although the theoretical understanding of the ground state is still poor, we 
do have very strong indirect evidence that Nambu's conjecture is right - we 
know why the energy gap of QCD is small. 

2 Lattice results for M n and F n 

As pointed out by Gell-Mann, Oakes and Renner [10] , the square of the pion 
mass is proportional to the strength of the symmetry breaking, 

Ml oc (m u + m d ) . 

This property can now be checked on the lattice, where - in principle - 
the quark masses can be varied at will. In view of the fact that in these 
calculations, the quarks are treated dynamically, the quality of the data is 
impressive. The masses are sufficiently light for %PT to allow a meaningful 
extrapolation to the quark mass values of physical interest. The results 
indicate that the ratio M% / (m u + ma) is nearly constant out to values of 
m u , m d that are about an order of magnitude larger than in nature. 

The Gell-Mann-Oakes- Renner relation corresponds to the leading term in 
the expansion in powers of the quark masses. At next-to-leading order, the 
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expansion in powers of m u , rrid (mass of the strange quark kept fixed at the 
physical value) contains a logarithm [11, 12]: 

M - A 4 + 3^ hl f +0(M4) }' W 

where M 2 = B(m u + m^) stands for the term linear in the quark masses. 
Chiral symmetry fixes the coefficient of the logarithm in terms of the pion 
decay constant F n , but does not determine the scale A 3 . A crude estimate 
for this scale was obtained more than 20 years ago [13], on the basis of the 
SU(3) mass formulae for the pseudoscalar octet. The result is indicated at 
the bottom of the left panel in figure 1. The other entries represent recent 
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Figure 1: Determinations of the effective coupling constants £3 and 



lattice results for this quantity [14]-[20]. The one of the RBC /UKQCD 
collaboration, £3 = 3.13 ± 0.33 stat ±0.24 syst [18], for instance, which concerns 
2+1 flavours and includes an estimate of the systematic errors, is considerably 
more accurate than our old estimate based on SU(3), £3 = 2.9 ± 2.4 [13]. 

The right panel shows the results for the scale A4, which determines the 
quark mass dependence of the pion decay constant at NLO of the chiral 
expansion. The analog of formula (1) reads 



F 1 



M 2 



16tt 2 F 2 



M 2 



(2) 



where F is the value of the pion decay constant in the limit m u ,md — > 0. 
A couple of years ago, we obtained a rather accurate result for A 4 , from 
a dispersive analysis of the scalar pion form factor [21]. The plot shows 
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that the lattice determinations of £4 have reached comparable accuracy and 
are consistent with the dispersive result. For a detailed discussion of the 
properties of the scalar pion form factor, I refer to [22]. This quantity is now 
also accessible to an evaluation on the lattice [23]. 



3 S-wave nn scattering lengths 

The hidden symmetry not only controls the size of the energy gap, but also 
determines the interaction of the Goldstone bosons at low energies, among 
themselves, as well as with other hadrons. In particular, as pointed out by 
Weinberg [24], the leading term in the chiral expansion of the S-wave tttt 
scattering lengths (tree level of the effective theory) is determined by the 
pion decay constant. The corresponding numerical values of a{] and are 
indicated by the leftmost dot in figure 2, while the other two show the result 
obtained at NLO and NNLO of the chiral expansion, respectively. The exotic 
scattering length is barely affected by the higher order corrections, but the 
shift seen in a® is quite substantial. The physics behind this enhancement 
of the perturbations generated by m u and is well understood: it is a 
consequence of the final state interaction, which is attractive in the S^-wave, 
rapidly grows with the energy and hence produces chiral logarithms with 
large coefficients. 




Figure 2: Comparing the theoretical predictions for the irn S-wave scattering 
lengths with lattice results (left) and with experiment (right) 



Near the center of the Mandelstam triangle, the contributions from higher 
orders of the chiral expansion are small [21]. Using dispersion theory to reach 
the physical region from there, we arrived at the remarkably sharp predictions 
for the two scattering lengths indicated on the left panel of figure 2. Our 
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analysis also shows that the corrections to Weinberg's low energy theorem 
for (2q, Qjq 3XG dominated by the effective coupling constants £ 3 , £ 4 discussed 
above - if these are known, the scattering lengths can be calculated within 
small uncertainties. Except for the horizontal band, which represents a direct 
determination of based on the volume dependence of the levels [25], all of 
the lattice results for the scattering lengths shown on the left panel of figure 
2 are obtained in this way from the corresponding results for £ 3 and £4. The 
figure demonstrates that the lattice results confirm the predictions for a[j, ag. 

4 Precision experiments at low energy 

The right panel of figure 2 compares the predictions for the scattering lengths 
with recent experimental results. While the K e 4 data of E865 [26], the 
DIRAC experiment [27] and the NA48 data on the cusp in K — ► 3n [28] 
all confirm the theoretical expectations, the most precise source of informa- 
tion, the beautiful K ei data of NA48 [29], gave rise to a puzzle. The Watson 
theorem implies that - if the electromagnetic interaction and the difference 
between m u and rrid are neglected - the relative phase of the form factors 
describing the decay K — > evnir coincides with the difference Sq — 5\ of scat- 
tering phase shifts. At the precision achieved, the data on the form factor 
phase did not agree with the theoretical prediction for the phase shifts. 



20 

15 

5°-5 l l10 
5 



-5 

0.28 0.3 0.32 0.34 0.36 0.38 0.4 

GeV 

Figure 3: Comparison of K e 4 data with the prediction for 5q — 8\ 

The origin of the discrepancy was identified by Colangelo, Gasser and 
Rusetsky [30]. The problem has to do with the fact that a K + may first 
decay into e + z/7r°7r° - the pair of neutral pions then undergoes scattering 
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and winds up as a charged pair. The mass difference between the charged 
and neutral pions affects this process in a pronounced manner: it pushes the 
form factor phase up by about half a degree - an isospin breaking effect, due 
almost exclusively to the electromagnetic interaction. 

Figure 3 shows that the discrepancy disappears if the NA48 data on 
the relative phase of the form factors are corrected for isospin breaking. 
Accordingly, the range of scattering lengths allowed by these data, shown 
on the right panel of figure 2, is in perfect agreement with the prediction. 
As indicated on the left panel of that figure, the low energy theorem for the 
scalar radius of the pion correlates the two S-wave scattering lengths to a 
narrow strip. If this correlation is used, the analysis of the K ei data leads to 
<2q = 0.220 ± 0.005 stat ± 0.002 syst [31]. The result has the same precision as 
the theoretical prediction and hits it on the nail. I conclude that the puzzle 
is gone: K e4: confirms the theory to remarkable precision. 

5 Expansion in powers of m s 

The examples discussed above all concern the effective theory based on 
SU(2)xSU(2), where the quantities of interest are expanded in powers of 
m u ,md, while m s is kept fixed at the physical value. The corresponding ef- 
fective coupling constants F, B,£i, . . . are independent of m u and ma, but 
do depend on m s . Their expansion in powers of m s can be worked out in 
the framework of the effective theory based on SU(3) xSU(3). For F, B, the 
expansion starts with 1 

F = F {l + ^i;-ji Jf + O(mJ)} , (3) 

B = B 1 1 + 1 -^(2L r 6 - LI) -fi v + 0{ml) } . 

The constants F , B represent the values of F, B in the limit m s — > 0. 
At NLO, only the coupling constants L 4) L Q of the chiral SU(3) xSU(3) La- 
grangian enter, weighted with the square of the kaon mass in the limit 
mu-> m d — 0, which I denote by M K . In this limit, the octet of Gold- 
stone bosons contains only three different mass values: M n = 0, Mx and 
M v . To the accuracy relevant in the above formulae, we have M\ = Bom s , 
= ^B m s . Up to corrections of higher order, M K and M v may be ex- 

1 The contributions of 0{m 2 s ) are also known explicitly, not only for F, B, but also for 
the coupling constants £i, . . . ,£7, which specify the effective Lagrangian at NLO [32, 33]. 
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pressed in terms of the physical masses as 



Ml = Ml-\Ml, Ml = \Ml-\Ml. (4) 

The chiral logarithms occurring in the above formulae may be expressed in 
terms of the function 

w- \ I- 

»* = P = K ^ < 5 > 

They involve the running scale \x at which the chiral perturbation series is 
renormalized, but the scale dependence of the renormalized coupling con- 
stants L4, Lg ensures that the expressions in the curly brackets of equation 
(3) are scale independent. 
The quark condensate, 

E=|(0|B«|0>| m ^ Q , (6) 

is determined by the same two constants: X = F 2 B. Accordingly, the ex- 
pansion of the condensate in powers of m s starts with 

S = S 1 1 + -y^Ll - 2n-- n n + 0(m 2 s ) j . (7) 

The coupling constants L 4 , L 6 as well as the loop graphs responsible for the 
chiral logarithms represent effects that violate the Okubo-Zweig-Iizuka rule. 
In the large N c limit, the quantities F,B,Y, become independent of m s , so 
that the ratios F/Fq, B/Bq, S/Xo tend to 1. If the OZI rule is a good guide 
in the present context, then these ratios should not differ much from 1. For 
a discussion of the implications of large OZI violations in these ratios, see 
[34]. The paramagnetic inequalities of Stern et al. [35] indicate that the sign 
of the deviations F/F — 1 and S/So — 1 is positive. 



6 Violations of the OZI rule ? 

Figure 4 compares recent lattice results for the dependence of the condensate 
on m s [14, 19, 36] with phenomenological estimates found in the literature 
[37, 38, 39, 40]. The latter are calculated from the values quoted for the 
running coupling constant L 6 , using the relation (7) with F = F n . The 
errors shown exclusively account for the quoted uncertainties in the coupling 
constants, while those arising from the corrections of 0(m 2 s ) are neglected 
(accounting for the difference between F and F n in the curly bracket of (7), 
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Figure 4: Violation of the OZI rule in the quark condensate 

for instance, generates corrections of this type). The plot shows that the 
uncertainties in the phenomenological estimates are large. Unfortunately, 
the lattice results are not yet conclusive, either. Note that some of these are 
preliminary and do not include an estimate of the systematic errors. 

In contrast to the lattice results for the condensate, those for the constant 
B are consistent with one another. As can be seen in figure 5, the values 
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Figure 5: OZI violations in the leading effective coupling constants 

obtained for B/Bq do not indicate a large violation of the OZI rule. This 
implies that the discrepancy seen in the lattice results for E = F 2 B originates 
in the factor F 2 . Indeed, the values quoted for F/F in [36] are puzzling, 
for the following reason. The quantity F„ represents the pion wave function 
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at the origin. The value of F K is somewhat larger, because one of the two 
valence quarks is heavier than in the case of the pion. Hence it moves mores 
slowly, so that the wave function is more narrow and thus higher at the 
origin: 2 

F K /F, = 1.192(7) . (8) 

If the value of F/F was larger than this, we would have to conclude that 
the wave function is more sensitive to the mass of the sea quarks than to the 
mass of the valence quarks. I do not see a way to rule this logical possibility 
out, but it is counter intuitive and hence puzzling. 

For the time being, the only conclusion to draw is that the lattice results 
confirm the paramagnetic inequalities and indicate that the constant B - the 
leading term in the expansion of ill 2 , in powers of m u and - does obey 
the OZI rule. Some of the data indicate that this rule approximately holds 
also in the case of F/F , but others suggest rather juicy violations in that 
case. The slow, but steady progress being made on the lattice gives rise to 
the hope that the dust will settle soon. 

7 Problems with scalar meson dominance ? 

The analysis of the lattice data on the quark mass dependence of M n , M K) F w , 
Fx in terms of the xPT representation to two loops will also make it possible 
to determine those coupling constants of the effective NNLO Lagrangian that 
contribute to these quantities. The theoretical estimates for those couplings 
[l]-[6], [44] rely on the assumption that the relevant sum rules are saturated 
by the lowest resonances. I know of two cases, where calculations within this 
framework run into a problem: 

• The data on nuclear (3 decay lead to a remarkably accurate value for the 
element V ud of the CKM matrix: V ud = 0.97418(26) [45]. The unitarity 
of this matrix then implies V us = 0.2258(11). Since the data on the rate 
of the semileptonic K decays require f+ * {0)V US = 0.21661(47) [46], 
the value of the form factor at the origin is determined very sharply: 
/«"°f-(0) = 0.9594(49). The recent lattice result, ° W J0) = 0.964(5) 
[47], as well as those for the ratio F K /F W [48], which offer an indepen- 
dent determination of V us , are consistent with this value. The chiral 

2 Note that, a few years ago, improved measurements performed at several laboratories 
demonstrated that the decay rate had been underestimated by more than 3 a. As 
a consequence, the values obtained for the CKM matrix element V us were too small, so 
that the ratio Fk/F^ of decay constants came out too large [41]. The value quoted in Eq. 
(8) is taken from Bernard and Passemar [42] . The mini- review of Rosner and Stone [43] , 
prepared for the 2008 edition of the PDG tables confirms this result within errors. 
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representation of the relevant form factors was worked out to NNLO by 
Bijnens and Talavera [49]. The evaluation of their representation for 
ff*~(0) with resonance estimates for the effective coupling constants 
[49, 50, 51] leads to values that are significantly higher than those 
above: the recent update of the calculation described in [52] leads to 
ff^(0) = 0.986 ± 0.007 1/JVc ± 0.002m 5 ,m p . 

• The branching ratios of the transitions K ± 7t ± tt , K° — > 7r + 7r~, 
K° — > 7r°7T are affected by the final state interaction. Conversely, the 
observed values of these ratios can be used to determine the phase 
difference between the two S- waves at s = M|-. In the past, work on 
this problem invariably led to a value for the phase difference that is too 
large, presumably because isospin breaking, which plays a crucial role 
here, was not properly accounted for. Only rather recently, Cirigliano, 
Ecker, Neufeld and Pich have performed a complete analysis of these 
transitions, based on %PT to NLO [53]. Their calculation accounts 
for isospin breaking, both from m u ^ and from the electromagnetic 
interaction [54]. Unfortunately, however, the discrepancy persists: their 
result is 5° (M 2 K ) - 5 2 {M 2 K ) = 57.5° ± 3.4° [55], while the Roy analysis 
of 7T7T scattering implies 3 5%(M 2 K ) - 5l{M 2 K ) = 47.7° ± 1.5° [21]. 

The above discussion of f^° n (0) assumes that, at the accuracy of interest, 
Ke 3 decay is properly described by the Standard Model, where the CKM 
matrix is unitary. Also, it relies on the value of V u d extracted from nuclear (3 
decay. These ingredients should not be taken for granted [56] , 4 but I consider 
it more likely that the discrepancy originates in the chiral calculation, also 
in the case of the phase difference 5q{M] < ) — 5q(M|-). I do not doubt the 
chiral representations of the form factors [49] and of the phase difference [53] , 
but the estimates used for the effective coupling constants play an equally 
important role. In my opinion, this is the weakest point in the above two 
applications of %PT . 

In the case of SU(2), the size of the coupling constants occurring in the 
NLO effective Lagrangian can be understood on the basis of vector meson 

comparison with more recent work on the 7T7t phase shifts is made in section 11. 
4 The results obtained from preliminary data on r decays into final states with strange- 
ness, for instance, do not agree with the above value of V us [57]. Also, a recent experiment 
on neutron decay [58] came up with a neutron lifetime that strongly disagrees with the 
world average. If confirmed, this calls for an increase in the value of gA/gv or in the value 
of V uc i- For the time being, the uncertainties in gA/gv are too large for neutron decay to 
compete with the superallowed nuclear transitions. In particular, the value of the neutron 
lifetime reported in [58] is consistent with the value of V u d obtained in [45] . 
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dominance alone 5 [13]. The above two examples, however, concern SU(3), 
where m s is also treated as a perturbation. In this case, the effects generated 
by the quark mass term in the Lagrangian of QCD are much more important. 
Since this term is a scalar operator, resonance estimates for those coupling 
constants that describe the dependence of the effective SU(3) Lagrangian on 
the quark masses rely on the scalar meson dominance hypothesis. In my 
opinion, it is questionable whether the complex low energy structure of the 
scalar states (strong continuum related to the rapidly rising ittt interaction, 
broad bump associated with the a, narrow peak from the /o(980), glueballs, 
etc.) can adequately be accounted for with this hypothesis. 

As shown in [60], the phenomenological results [37] for the coupling con- 
stants occurring in the SU(3) Lagrangian at NLO can be understood on the 
basis of the assumption that the contributions from the lowest resonances 
with spin < 1 dominate. At first sight, this might appear to confirm the va- 
lidity of scalar meson dominance in the present context, but that is not the 
case: in [60] , the quark mass dependence of the effective Lagrangian is taken 
from phenomenology (the model used for the scalar resonances contains two 
free parameters and these are tuned so as to reproduce the observed values 
of L 5 and L 8 ). Phenomenological information about the quark mass depen- 
dence of the NLO couplings relevant for K — > nrr is lacking. In this case, 
factorization is used to estimate the effective coupling constants [61]. The 
contributions generated by a single resonance do factorize, but in the more 
complex situation encountered in the scalar channel, factorization may fail. 

In the first example, the available experimental and theoretical informa- 
tion should suffice to determine all of the coupling constants relevant for the 
form factors to NNLO. Those accounting for the slope and the curvature 
were worked out already [50]. It would be very instructive to determine the 
remaining ones, which describe the dependence of the form factors on the 
quark masses. Inserting the values obtained in this way in the chiral repre- 
sentation of the form factors should lead to a coherent picture, but the results 
obtained for some of the coupling constants will necessarily differ from the 
resonance estimates. It would be most useful to understand the origin of 
the difference, because similar departures from resonance saturation must be 
expected in other matrix elements. Whether this will lead to a resolution of 
the second discrepancy remains to be seen - at any rate, the puzzle with the 
phase determination via K — > ixix cries for a resolution. 

5 The more sophisticated treatment of the vector meson dominance hypothesis described 
in [59] confirms the crude estimate for £4 given in equations (C.12), (C.13) of [13]. 
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8 Puzzling results in decay 



The low energy theorem of Callan and Treiman [62] predicts the size of the 
scalar form factor of the decay K — > nfiu at one particular value of the 
momentum transfer, namely t = M\ — M 2 : 

f (M 2 K - Ml) = ^ + 0(m u , m d ) . (9) 

Within QCD, the relation becomes exact if the quark masses m u and md 
are set equal to zero. The corrections of first nonleading order, which have 
been evaluated long ago [63], are tiny: they lower the right hand side by 
3.5 x 10~ 3 . In the meantime, the chiral perturbation series of fo{t) has been 
worked out to NNLO [49, 64]. As pointed out by Jamin, Oiler and Pich [50], 
the curvature of the form factor can be calculated with dispersion theory. 
Their dispersive representation agrees very well with the more recent one 
of Bernard, Oertel, Passemar and Stern [65]: theory reliably determines the 
curvature of the form factor. Accordingly, the theoretical prediction for the 
value at the Callan- Treiman point, t = M\ — can be converted into 
a prediction for the slope. The result obtained by Jamin et al. in 2004 was 
Ao = 0.016(1). The update of their calculation with the improved information 
available in 2006 led to A = 0.0147(4). Within the remarkably small errors, 
this agrees with the outcome of the recent analysis described in [42] , for which 
the Standard Model prediction is Ao = 0.0150(8) [66]. 

Recently, the NA48 collaboration published their results for the K® 3 form 
factors [67]. Their result for the scalar slope, A = 0.0117(7)(1), is in flat 
contradiction with the theoretical prediction just discussed. 

The NA48 experiment is not the first to measure the slope of the scalar 
form factor relevant for decay. Figure 7 compares the outcome of this 
experiment with results obtained by ISTRA [68], KTeV [69], KLOE [70] 
and with earlier findings, taken from the PDG listings of 2004 [71]. The 
earliest result shown indicates the outcome of the high statistics experiment 
of Donaldson et al. [72], which came up with a slope of A = 0.019(4) and 
thereby confirmed the validity of the Callan- Treiman relation. The plot 
shows that quite a few of the experiments performed since then obtained 
quite different results. 

In 1985, when we worked out the corrections to the Callan- Treiman rela- 
tion to one loop of %PT [63] , the experimental situation was entirely unclear. 
We emphasized that there is no way to reconcile some of the published ex- 
perimental results with the Standard Model. In the meantime, both the 
accuracy of the theoretical prediction and the quality of the data on the 
scalar form factor improved considerably, but an experimental discrepancy 



13 



1.3 
1.25 
12 

f Q (t) 1,5 

1.1 
1.05 
1.0 

0.05 0.1 0.15 0.2 j 0.25 

t (GeV 2 ) M ^ - M* 

Figure 6: NA48 results for the scalar K^ 3 form factor 

persists: while the results of ISTRA, KTeV and KLOE are consistent with 
one another, they disagree with NA48. 

There are not many places where the Standard Model fails at low ener- 
gies. Hints at such failures deserve particular attention. As pointed out by 
Jan Stern and collaborators [65], the Callan-Treiman relation can serve as 
a probe for the presence of right-handed quark couplings of the VT-boson - 
the existing limits on those couplings are not very stringent. It is premature, 
however, to interpret the NA48 data as evidence for the occurrence of effects 
beyond the Standard Model. For the time being, we are merely faced with 
an experimental discrepancy. The experiment is difficult, because the tran- 
sition is dominated by the contribution from the vector form factor, while 
the Callan Treiman relation concerns the scalar form factor. In particular, 
the radiative corrections must properly be accounted for [52, 73, 74, 75] . As 
emphasized by Franzini [76], the data analysis must cope with very strong 
correlations between the slopes of the vector and scalar form factors. Also, 
for the analysis of the data to match their quality, it is essential that the con- 
straints imposed by dispersion theory are respected - publishing fits based 
on linear parametrizations of the t-dependence or on pole models is meaning- 
less. An analysis of the charged kaon decays collected by NA48 might help 
clarifying the situation. 
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Figure 7: Old and new results for the slope of the scalar form factor 

9 Dispersion theory 

In the remainder of these lecture notes, I discuss some of the progress made 
in carrying the low energy analysis to higher energies. First steps in this 
direction were taken by Gasser and Meissner [77], who compared the rep- 
resentation of the scalar and vector form factors obtained in the framework 
of xPT to two loops with the dispersive representation. In particular, they 
determined the range of validity of the representation obtained by truncating 
the chiral perturbation series at one or two loops. Indeed, many of the issues 
discussed in the first part of these lectures involve dispersion theory: the con- 
straints imposed on the form factors by analyticity and unitarity (Watson 
final state interaction theorem) play a central role in the data analysis. 

In the following, I consider a different example: the tttt scattering ampli- 
tude. The chiral perturbation series of this amplitude, even if truncated only 
at NNLO, is useful only in a very limited range of the kinematic variables 
- definitely, the poles generated by p-exchange are outside this region. The 
range can be extended considerably by means of dispersion theory, exploit- 
ing the fact that analyticity, unitarity and crossing symmetry very strongly 
constrain the low energy properties of the scattering amplitude. 

From the point of view of dispersion theory, scattering is particularly 
simple: the s-, t- and li-channels represent the same physical process. As 
a consequence, the scattering amplitude can be represented as the sum of 
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a subtraction term and a dispersion integral over the imaginary part. The 
subtraction term involves two subtraction constants, which may be identified 
with the two S-wave scattering lengths. The dispersion integral exclusively 
extends over the physical region [78]. 

The projection of the amplitude on the partial waves leads to a dispersive 
representation for these, the Roy equations. I denote the S-matrix elements 
by Si = T)\ exp 2ib\ and use the standard normalization for the corresponding 
partial wave amplitudes t\: 



S I e (s) = l + 2ip(s)t I i (s), p(s) = y/l-4M*/s. (10) 

The S-matrix elements and the partial wave amplitudes are analytic in the 
cut s-plane. There is a right hand cut (4M 2 < s < oo) as well as a left hand 
cut (—00 < s < 0). The Roy equation for the partial wave amplitude with 
/ = £ — 0, for instance, reads 

" ^ poo 

f°(s)=a + ( S -4ilO& + ££ / ds'K^^lmt'is'). (11) 

7=0 i=0 JiM i 



As mentioned above, the equation contains two subtraction constants, which 
can be expressed in terms of the S-wave scattering lengths: 



b = (2a° - 5a 2 )/12Ml (12) 



The kernels K\f, (s, s') are explicitly known algebraic expressions which only 
involve the variables s, s' and the mass of the pion, e.g. 

00 1 2 hvjjs + sf- AMl)/s>} 5s' + 2s- 

001 S,S) ~ jr(s'-s\ + ^ 3tt( S -4M, 2 ) ^ 3tt£(£ - 4M, 2 ) / [ 6) 

r.h.cut l.h.cut 

The integrals on the right hand side of (11) thus only involve observable 
quantities: the imaginary parts of the partial waves. 

The pioneering work on the physics of the Roy equations was carried out 
more than 30 years ago [79]. The main problem encountered at that time 
was that the two subtraction constants were not known. These dominate the 
dispersive representation at low energies, but since the data available at the 
time were consistent with a very broad range of S-wave scattering lengths, the 
Roy equation analysis was not conclusive. The insights gained by means of 
%PT thoroughly changed the situation. Since the S-wave scattering lengths 
are now known very accurately, the Roy equations have become a very sharp 
tool for the analysis of the nrr scattering amplitude. 
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10 Mathematics of the Roy equations 



The mathematical properties of the Roy equations are quite remarkable and 
are discussed in detail in the literature. For an extensive review and refer- 
ences to the original literature, I refer to [80]. In the following, I restrict 
myself to those features that are essential for an understanding of the con- 
sequences of these equations. It is convenient to indicate the isospin of the 
partial waves as an index: S°, S 2 denote the S-waves of isospin and 2, 
respectively, the P-wave is referred to as P 1 , while D°, D 2 denote the two 
D- waves etc. 

The Roy equations were derived from axiomatic field theory, for values 
of s in the interval — 4M 2 < s < 60M 2 [78]. If the scattering amplitude 
obeys Mandelstam analyticity, then the derivation goes through on a slightly 
larger domain, — 4M 2 < s < 68M 2 , but even then, we can make use of these 
equations only on a finite interval. The upper end of the interval on which 
the Roy equations are solved is referred to as the matching point. I denote 
the corresponding value of s by s m . 

If we for the moment treat the imaginary parts of all other partial waves 
and the two scattering lengths as a given input, the Roy equation (11) 
amounts to a representation of the function Retp(s) as an integral over 
Imtg(s') plus a known remainder. In the elastic region, unitarity imposes 
a second such relation - we thus have two equations for the two unknown 
functions Ret[j(s) an d Imt[j(s). Hence no freedom appears to be left: book- 
keeping suggests that equation (11) determines these two functions in terms 
of the given input. 

If we restrict ourselves to the elastic region (s m < 16M 2 ), then the naive 
expectation is indeed correct: the solution is unique. In order to push the 
matching point to higher energies, we need to know the elasticity ?7q(s), but if 
that is the case, we still have two equations for two unknowns. The properties 
of the system of equations, however, change if the matching point is pushed 
up. The solution remains unique only if the phase at the matching point stays 
below jir. If #o(s m ) * s m the interval \n < 5o(s m ) < vr, the system admits 
a one-parameter family of solutions, for n < <5g(s m ) < |"7r, the manifold of 
solutions is of dimension 2, etc. The available phase shift analyses indicate 
that the phase 5q(s) passes through \-k somewhere around 0.85 GeV, goes 
through 7r in the vicinity of KK threshold, but stays below §7T on the entire 
range where the Roy equations are valid. If the matching point is taken 
at the upper end of this range, equation (11) thus admits a two-parameter 
family of solutions [81]. 

Analogous statements also hold for the Roy equations obeyed by the other 
partial waves. For the P-wave, the phase 5\{s) reaches |7r around the mass 
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of the p and remains below n for s < 68M^. Accordingly, if the matching 
point of the P-wave is taken above M p , the Roy equation for t\(s) admits 
a one-parameter family of solutions. The exotic S-wave, Sq(s), on the other 
hand, is negative on the entire range where the Roy equations are valid and 
stays above In this case, the number of free parameters is equal to -1, 

irrespective of the choice of the matching point: the Roy equation for £q(s) 
does not in general admit a solution if the imaginary parts of all other waves 
(as well as Imt^s) for s > s m ) are prescribed arbitrarily. The input must 
be tuned in order for a solution to exist at all. For all of the higher partial 
waves, the phase remains below |7T on the interval where the Roy equations 
hold, so that the solution is unique - if the phase at the matching point is 
negative, the input needs to be tuned for a solution to exist. 

The higher waves only play a minor role at low energies. Although, in 
principle, their properties are correlated with those of the S- and P-waves, 
we may first use a phenomenological parametrization for the partial waves 
with £ > 2 and solve the Roy equations for the S- and P-waves with this 
input. Then, the Roy equation for the D-, F-, G-waves can be solved, one 
by one, using the representation for the S- and P-waves obtained in the first 
step (on the range where the Roy equations are valid, the contributions from 
the imaginary parts of the partial waves with £ > 5 are too small to matter). 
If the S- and P-waves are known, the Roy equations for the higher waves 
fix their behaviour at low energies within very small uncertainties. Finally, 
we may perform an iteration, inserting the representation obtained for the 
imaginary parts of the waves with £ > 2 in the Roy equations for the S-and 
P-waves. There is no need for further iterations because the changes found 
in these waves are tiny. 

The heart of the matter is a coupled system of three integral equations 
for the three partial waves £q(s), t\(s), tl(s). The input of the calculation 
consists of the following parts: 

1. S-wave scattering lengths 

2. Elasticities of the S-and P-waves below s rn 

3. Imaginary parts of the S- and P-waves above s rn 

4. Imaginary parts of the higher partial waves 

If the matching point is taken below the p mass, the system only admits a 
solution if this input is tuned. If the matching point is somewhere in the range 
between 0.77 and 0.85 GeV (more precisely, the range where 5{(s m ) > |7r, 
but 5g(s m ) < \ti), the solution is unique. Above, 0.85 GeV, but below 2Mk, 
there is a one-parameter family of solutions and for s rn > 4M|-, in particular, 
if the Roy equations are solved on the entire interval where they are valid, 
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the solution contains two free parameters [82] . In order to arrive at a unique 
solution, we may, for instance prescribe two of the phases at a suitable energy. 
It is convenient to use the values of Sq(s) and 5{(s) at y^s = 0.8 GeV for this 
purpose. I denote the value of s at this energy by 

s A = (0.8 GeV) 2 . 

11 Low energy analysis of 7T7T scattering 
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Figure 8: Low energy behaviour of the partial wave amplitude i°(s) 

At low energies, the two S-wave scattering lengths are the most important 
parameters, because there, the right hand side of the Roy equations is domi- 
nated by the subtractions. Figure 8 demonstrates that the subtraction term 
in equation (11) dominates the behaviour of the partial wave t§ (s) through- 
out the region shown, which extends to about 500 MeV. As predicted by 
current algebra (tree graphs of the effective theory), ig(s) contains an Adler 
zero and grows approximately linearly with s. 

As discussed in the preceding section, the behaviour in the low energy re- 
gion is not fully determined by the input listed, but involves two additional 
low energy degrees of freedom, which can be identified with 5q(sa), SI(sa)- 
The value of 5{(s A ) is known well, because the Watson final state interaction 
theorem connects the P-wave phase shift with the vector form factor of the 
pion, which is accurately determined by the beautiful data on e + e~ — > 7r + 7r~ 
and t~ -> v T The result obtained in [80] reads 5\(s A ) = 108.9° ± 2°. 
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The experimental information about 8q(sa), on the other hand, is compar- 
atively meagre - this currently represents the main source of uncertainty in 
low energy 7nr scattering. In [80], we observed that phase differences are more 
easily determined than the phases themselves. Indeed, in contrast to the ex- 
perimental information about 5q(sa), the one for S\(sa) — Sq(sa) * s perfectly 
coherent. Together with the value of SI(sa) quoted above, the experimental 
information about this difference implies Sq(sa) = 82.3° ± 3.4°. 

In order to explore the space of solutions, we first fixed the matching 
point at 0.8 GeV and ignored the predictions of xPT, treating free 
parameters [80]. The consequences of the low energy theorems for Oq,Oq 
were worked out separately, in [21], where it was shown, for instance, that 
the scattering lengths and the effective ranges of the S-, P-, D- and F- waves 
can be calculated within remarkably small uncertainties on this basis. 

One of the sources of uncertainty arises from the "high energy" part 
of the input: the dispersion integrals extend to infinity - phenomenology is 
used to estimate the contributions from the region above the matching point. 
For the asymptotic properties of the amplitude, we relied on the literature, 
in particular on the work of Pennington [83], who had carefully examined 
the relationship between the behaviour at low energies and the high energy 
properties of the amplitude. For an update of the Regge parametrization 
used in this context, see [84]. It is essential that the Roy equations involve 
two subtractions, so that the kernels fall off with the third power of the 
variable of integration. The left hand cut plays an important role here: as 
can be seen in equation (13), the part of the kernel that accounts for the right 
hand cut falls off only with the first power of the variable of integration, but 
the high energy tail is cancelled by the contribution from the left hand cut. 
This ensures that the contributions from the low energy region dominate. 

The results obtained in [80] were confirmed by Stern and collaborators 
[85]. Moreover, these authors made fits to the data available at the time, in 
order to determine the scattering lengths from experiment and to compare the 
results with the theoretical predictions. The experimental information about 
the low energy behaviour of S 2 played an important role in their analysis. 
Unfortunately, this information is not coherent. In particular, the authors 
had to make a choice between the two phase shift analyses reported in [86]. 
With the choice made, the results obtained for a[j, a 2 , turned out not to be 
in satisfactory agreement with the theoretical predictions. A similar analysis 
has now been carried out for the K e4 data of NA48/2 [29, 31] - as mentioned 
already in section 4, these data confirm the theory to remarkable precision. 

The approach of Yndurain and collaborators [87, 88] is quite different. 
These authors do not make an attempt at solving the Roy equations, but 
only use these to test and improve their parametrization of the data, making 
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fits where the difference between the left- and right hand sides of the Roy 
equations and the one between the parametrization and the data is given 
comparable weight. Indeed, their parametrization underwent a sequence of 
gradual improvements, which removed some of the deficiencies pointed out 
in [89, 90, 91]. In particular, the residue of the p-trajectory increased step 
by step, as a consequence of the Olsson sum rule. The representation of 
the higher partial waves was improved and the experimental determination 
of the scattering lengths described in KPY III [88] is now consistent with 
our predictions. Unfortunately, however, an important difference to our rep- 
resentation persists: their phase 5q(s) still contains a kink (discontinuity 
in the first derivative) at 932 MeV, as well as a hump below that energy. 
As discussed in detail in [90], these phenomena are artefacts, produced by a 
parametrization that is not flexible enough. The analysis in [92] corroborates 
this conclusion. 

Incidentally, the kink and the hump are also responsible for the remaining 
disagreement in some of the threshold parameters: the contributions from 
Imt[j(s) to the sum rules for these quantities are not the same (the sum rules 
are listed in Eqs. (14.1) and (14.3) of [80]). Replacing the parametrization 
of 5q(s) below KK threshold by ours, but taking all other contributions 
from KPY III [88], the result reproduces all of the entries listed in table 2 
of [21], within errors. In other words, as far as the integrals relevant for 
the threshold parameters are concerned, the only remaining difference that 
matters concerns the behaviour of 5q(s) below KK threshold. Work aimed 
at improving the quality of that part of the representation is in progress [93] . 

12 Behaviour of the S-wave with 1 = 

As mentioned in the preceding section, the value of Sq(sa) currently repre- 
sents the main source of uncertainty in low energy nn scattering. The quoted 
range, Sq(sa) = 82.3° ± 3.4°, which follows from the experimental informa- 
tion on the phase difference 5{(sa) — 5q{sa) an d on 8{(sa), does not cover 
all of the data on Sq(sa), which are contradictory. While the result of the 
experiment described in the 1973 PhD thesis of Wolfgang Ochs [94, 95], for 
instance, is perfectly consistent with the Roy equations and does lead to a 
value of Sq(sa) in this range, the phase shift analysis of the polarized data 
of the CERN-Cracow-Munich collaboration [96] calls for a higher value. As 
shown by Kamihski, Lesniak and Loiseau [97], the phase ambiguity occurring 
in that analysis can be resolved by means of the Roy equations. The solu- 
tions obtained by these authors are of good quality: the difference between 
input and output for the real parts are of order 10~ 3 , for the S-waves as 
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well as for the P-wave. According to figure 2a in [97], the resulting fit yields 
5q(sa) — 87°. In view of the relatively large errors attached to the phase shift 
in [96], this result must come with a sizable uncertainty and may thus not 
be inconsistent with the range obtained in [80], but it is on the high side. 
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Figure 9: Behaviour of Sq below KK threshold 

The parametrizations of Kamihski, Pelaez and Yndurain [88] yield even 
higher values: 5°(s A ) = 90.7° ± 0.7 (A), 5°(s A ) = 90.5° ± 0.7 (B). In view of 
the remarkably small error, these results disagree with those obtained from 
6{(sa) — ^o(sa) [80] or from a Roy equation fit to the data of [94]. One of 
the reasons for arriving at such a high value is that the authors include the 
result for the phase difference <5q(M|-) — $j(Mjj-) obtained from K — > tttt [53] 
in their fitting procedure. This pulls the value of 5q(s a ) up. The response of 
the Roy equations to this change in the input value for S^(sa) is an increase in 
8^{M 2 K )-5l{M 2 K ) of 2°. The fit obtained in KPYIII yields a somewhat larger 
shift: the value for <^(Mf) - <5o(M£) is 50.9° ± 1.2°, higher than our result 
by 3.2°. The difference is produced by the kink mentioned in the preceding 
section, which can also be seen in figure 9. The kink generates a violation 
of causality and hence of the Roy equations: while our amplitude or the 
one of Kamihski, Lesniak and Loiseau [97] do represent decent approximate 
solutions of the Roy equations, the one in KPYIII does not: in the region 
between 0.7 and 1 GeV, the difference between input and output for the real 
parts of the S- waves is of order 0.1. Quite irrespective of these details, the 
increase in the phase difference 8q(Mk) — ^(M^) produced by an increase 
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in the value of Sq(sa), even a very large one, is much too small to bridge 
the gap between scattering and decay, which is of the order of 10°. The 
puzzle discussed in section 7 needs to be resolved before information about 
mi scattering can reliably be extracted from the decay K — > ixix . 

Although the possibility that the phase Sq(sa) is above the range obtained 
from S\(sa) — Sq(sa) looks unlikely, it cannot be ruled out entirely. For this 
reason, when analyzing the pole position of the a [98], we stretched the error 
bar towards higher values and used 

6 (s A ) = 82.3^° . (14) 

The band in figure 9 shows the corresponding range of Roy solutions, together 
with the parametrizations of Sq(s) proposed in [88, 99, 100, 101]. The plot 
shows that the range (14) covers the values of 5q(sa) obtained with these. 

13 Pole formula 

The positions of the S-matrix poles represent universal properties of QCD, 
which are unambiguous even if the width of the resonance turns out to be 
large, but they concern the non-perturbative domain, where an analysis in 
terms of the local degrees of freedom - quarks and gluons - is not in sight. 
One of the reasons why the values for the pole position of the / (600) = a 
quoted by the Particle Data Group cover a very broad range is that all 
but one of these either rely on models or on the extrapolation of simple 
parametrizations: the data are represented in terms of suitable functions 
on the real axis and the position of the pole is determined by continuing 
this representation into the complex plane. If the width of the resonance 
is small, the ambiguities inherent in the choice of the parametrization do 
not significantly affect the result, but the width of the a is not small. For 
a thorough discussion of the sensitivity of the pole position to the freedom 
inherent in the choice of the parametrization, I refer to [92]. 

The determination of the a pole provides a good illustration for the 
strength of the dispersive method and for the relative importance of the 
various terms on the right hand side of the Roy equations. Using known re- 
sults of general quantum field theory [102, 103], we have shown that the Roy 
equations also hold for complex values of s, in the intersection of the relevant 
Lehmann-Martin ellipses [98]. The pole sits on the second sheet, which is 
reached from the first by analytic continuation from the upper half plane into 
the lower half plane, crossing the real axis in the interval 4M^ < s < 16M%, 
where the scattering is elastic. For a real value of s on this interval, we have 
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Figure 10: Domain of validity of the Roy equations (s in units of Ml) 

S I e (s±ie) 1 = exp±2i#(s), so that S^s+ie) 11 = Sj(s-ie) 1 = 1/ S^s+ie) 1 . 
Hence the relation 

Sl{s) H = l/SM 1 (15) 

holds on a finite interval of the real axis. Since the equation connects two 
meromorphic functions, it also holds for complex values of s. A pole on the 
second sheet thus occurs if and only if S\ (s) has a zero on the first sheet. The 
net result of this discussion is that we have an exact formula for resonances: 

5/(a)=0. (16) 

For resonances with the quantum numbers of the vacuum, I — £ — 0, the 
element Sq(s) of the S-matrix is relevant. It is specified explicitly in equations 
(10) and (11), in terms of the imaginary parts of the partial waves on the real 
axis. The formula thus exclusively involves observable quantities and can be 
evaluated for complex values of s just as well as for real values. It provides 
for the analytic continuation necessary to determine the resonance position 
- a parametrization is not needed for this purpose. 

14 The lowest resonance of QCD 

Inserting our central representation for the scattering amplitude in (11), we 
find that, in the region where the Roy equations are valid, the function Sq(s) 
has two zeros in the lower half of the first sheet: one at y/s = 441 — i 272 
MeV, the other in the vicinity of 1 GeV [98]. While the first corresponds to 
the a, the second zero represents the well-established resonance / (980). Our 
analysis sheds little light on the properties of the latter, because the location 
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of the zero is sensitive to the input used for the elasticity t)q(s) - the shape of 
the dip in t)q(s) and the position of the zero represent two sides of the same 
coin. For this reason, I only discuss the a. 

We are by no means the first to find a resonance in the vicinity of the 
above position. In the list of papers quoted by the Particle Data Group [104], 
the earliest one with a pole in this ball park appeared more than 20 years 
ago [105]. What is new is that we can perform a controlled error calculation, 
because our method is free of the systematic theoretical errors inherent in 
models and parametrizations. For this purpose, it is convenient to split the 
right hand side of the Roy equation for t®(s) into three parts: 

a. Subtraction constants 

b. Contribution from Imt^s) below KK threshold 

c. Contributions from higher energies and other partial waves 

a. Subtraction constants 

The subtraction constants are determined by the S-wave scattering lengths. 
The predictions for these read a° = 0.220 ±0.005 and a 2 = -0.0444 ± 0.0010 
[21]. Following error propagation, we find that an increase in a[j by 0.005 
shifts the pole position by (—2.4 + 2 3.8) MeV, while the response to an in- 
crease in dp by 0.0010 is a shift of (0.8 — i 4.0) MeV [98]. These numbers show 
that the error in the pole position due to the uncertainties in the subtraction 
constants are small. 

b. Contribution from ImtJ(s) below KK threshold 

Below KK threshold, the S-waves are elastic to a very good approxi- 
mation. As shown in figure 11, the function Imt[j(s) shows a broad bump, 
nearly hits the unitarity limit somewhere between 800 and 900 MeV and then 
rapidly drops, because the phase steeply rises, reaching 180° in the vicinity 
of 2 Mk- Hence there is a pronounced dip in Imt^s) near KK threshold. 
The bump seen in the imaginary part in figure 11 may be viewed as a picture 
of the broad resonance we are discussing here. 

Below KK threshold, the behaviour of the imaginary part is controlled 
almost entirely by the phase shift <5q(s). Replacing the integral over our 
central representation for Imt[](s) from 4M^ to AM"^ by the one of Bugg 
[99], but leaving everything else as it is, the pole moves to 444 - i 267 MeV. 
Repeating the exercise with the representations of Kamihski, Pelaez and 
Yndurain [88], Achasov and Kiselev [100] and Albaladejo and Oiler [101], 
the pole is shifted to 458 - i 253 MeV, 438 - i 274 MeV and 451 - i 257 MeV, 
respectively. If the imaginary part of t[J(s) is evaluated with the lower edge 
of the band shown in figure 9, which corresponds to £q(sa) = 78.3°, the pole 
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Figure 11: Behaviour of Imt^ below KK threshold 

occurs at 435 - i 276 MeV, while for the upper edge of the band, characterized 
by 5%(s A ) = 92.3°, the pole sits at 456 - i 262 MeV. 

c. Contributions from higher energies and other partial waves 

Finally, I turn to the contributions of the third category: higher energies 
and other partial waves. Among these, the one from the P-wave, for example, 
is by no means negligible, but, as mentioned above, this wave is known very 
well. In fact, in the vicinity of the zero of Sq(s), the sum of the contributions 
of this category can be worked out quite accurately. The net uncertainty 
in the pole position from this source is ±4 ± i 6 MeV. As a check, we can 
simply replace our central representation for the contributions of category c 
by the one in [88], retaining our own representation only for the remainder. 
The operation shifts the pole position by —0.6 — i 1.2 MeV, well within the 
estimated range. 

Adding the errors up in square, the result for the pole position becomes 

= 441 +s 6 - i 272 + J l2 .5 MeV [98]. (17) 

The error bars account for all sources of uncertainty and are an order of 
magnitude smaller than for the estimate yfs^ = (400 - 1200) - i (250 - 500) 
MeV quoted by the Particle Data Group [104]: the position of the lightest 
resonance of QCD can now be calculated reliably and quite accurately. 
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